Oscillating tails of dispersion-managed soliton 
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Oscillating tails of dispersion-managed optical fiber system are studied for strong dispersion map in 
the framework of path-averaged Gabitov-Turitsyn equation. The small parameter of the analytical 
theory is the inverse time. An exponential decay in time of soliton tails envelope is consistent 
with nonlocal nonlinearity of Gabitov-Turitsyn equation, and the fast oscillations are described by a 
quadratic law. The pre-exponential modification factor is the linear function of time for zero average 
dispersion and cubic function for nonzero average dispersion. 



INTRODUCTION 

A dispersion management is key current technology for 
development of ultrafast high-bit-rate optical communi- 
cation lines H s n s nB s s s [11 m E m. 

A dispersion-managed (DM) optical fiber is designed to 
achieve low (or even zero) path-averaged group velocity 
dispersion by periodic alternation of the sign of the dis- 
persion along an optical line which dramatically reduces 
pulse broadening. 

If nonlinearity is neglected it would be possible to 
achieve transmission of optical pulses through DM sys- 
tem without significant distortion. However, the lin- 
ear transmission is limited by the nonlinear distance 
z < Zni = (cr|MoP)^^, which is determined by the Kerr 
nonlinearity a and the characteristic pulse power |uoP; 
z is the propagation distance along fiber. The charac- 
teristic power can not be chosen too small in order to 
maintain appropriate value of the signal-to-noise ratio. 
As a result, Kerr nonlinearity is essential for pulse prop- 
agation which results in distortion of optical pulses and 
finite bit error rate (BER) for information transmission. 
There are two alternative approaches to deal with optical 
nonlinearity. 

First approach is to reduce nonlinearity i.e. to make 
fiber system as close as possible to linear system. E.g. 
one can use low nonlinearity fiber [l^ Il5l| with increased 
cross section. Another approach is to compensate nonlin- 
earity (nonlinearity management) using either semicon- 
ductor devices |l6| or interferometer-based compensa- 
tion of nonlinear phase shift 0|. These systems usually 
refer to as quasi-linear systems and they attracted very 
much attention in recent years. Disadvantage of quasi- 
linear systems is that nonlinearity is still essential at long 
enough distance^ 0| (e.g. transoceanic distance) and 
it results in information losses. Also quasi-linear sys- 
tems often require complex and cumbersome dispersion 
compensation [l^l . 

Second approach is not to reduce nonlinearity but 
rather to use nonlineariy to achieve high bit rate trans- 
mission. This approach is based on a concept of using 
of DM solitons as the carrier of a bit of information. 



DM soliton is a solitary pulse which propagates in optical 
fiber [4, 5, 6, 7j.8., DM soliton experiences pe- 

riodic oscillations, as a function of propagation distance, 
with period of dispersion compensation. Single DM soli- 
ton can propagate in optical fiber without distortion (ex- 
cept periodic oscillations) for arbitrary long distance even 
for nonzero path-averaged dispersion (dispersion, aver- 
aged over period of dispersion variation). Propagation 
without distortion results from balance between nonlin- 
earity and dispersion on the period of dispersion compen- 
sation. In contrast, quasi-linear pulse would increase its 
width with propagation distance (in addition to periodic 
oscillations) due to nonzero path-averaged dispersion. 

DM soliton-based transmission has its own disadvan- 
tages caused by interaction between different DM soli- 
tons (representing different bits of information). There 
are two types of interactions. First type is caused by the 
interaction between DM solitons with the same carrier 
frequency (the same information channel). In addition, 
modern optical lines use a wavelength-division multi- 
plexing (WDM) which allows the simultaneous transmis- 
sion of several information channels, modulated at differ- 
ent wavelengths, through the same optical fiber. WDM 
means that there is also a second type of interaction, 
which is the interaction between DM solitons in different 
channels. Traditionally, the second type is considered to 
be the most dangerous because DM solitons in different 
channels move with different group velocity and, respec- 
tively, they pass through each other causing collisions 
through nonlinear interaction. Collisions result in a jit- 
ter in DM soliton arrival times and finite BER. However, 
recently, the record-breaking, 1.09Tbits/s, transmission 
was achieved at 18,000 km optical line based on DM soli- 
tons ^3 which brought increasing interest to DM soliton- 
based transmission systems. The interchannel interac- 
tions in that transmission experiment were suppressed 
by a periodic-group-delay dispersion-compensating mod- 
ule [l9|, and intrachannel interaction is most important. 
The main purpose of this Article is to determine the tails 
DM soliton which determines the interaction of DM soli- 
tons in the same channel. Soliton tails are responsible for 
interaction of pulses launched into optical fiber which es- 



sentially limits bit-rate capacity of optical line and makes 
finding DM soliton asymptotic behaviour an important 
practical and fundamental problem. It is found that 
envelope of soliton decays exponentially in time (much 
slower decay compare to Gaussian) so one can expect a 
strong interaction of different DM solitons. 



PATH-AVERAGED EQUATION 

Neglecting polarization effects, higher order dispersion, 
losses, stimulated Raman scattering and Brillouin scat- 
tering, the propagation of optical pulse in a DM fiber 
is described by a scalar nonlinear Schrodinger equation 
(NLS) 

iAz~]^p2{Z)ATT + ^{Z)\A\'A^Q, (1) 

where A is the envelope of optical pulse, Z is the prop- 
agation distance, T is the retarded time, P2{Z) is the 
group- velocity dispersion which is the periodical function 
of Z, (t{Z) ~ {2nn2)/ (^XqAej f{z)) is the nonlinear coeffi- 
cient, n2 is the nonlinear refractive index, Aq = 1.55 /im 
is the carrier wavelength, Aefj is the effective fiber area 
that in general case depends on Z. Typical values of dis- 
persion and nonlinear coefficient are (32 — — 20.0ps^/km, 
a = 0.0013 (kmmW)^^ for standard monomode fiber, 
and typical values for dispersion compensating fiber are 
(32 = 103.9 psVkm, a = 0.00405 (kmmW)"!. 

It is convenient to introduce the dimensionless vari- 
ables: z — Z/lo, t = T/to, u ~ A/\/Pq, where is the 
typical period of dispersion variation, is the typical 
pulse duration and Pq is the typical pulse power. For 
typical optical lines /o ^ 50 km, ^ lOps, Pq ^ 1 mW. 
Eq. in dimensionless variables takes the following 
form: 



iuz + d[z)utt + c(z)|m|^u — 0, 



(2) 



where c{z) = cr(z)Po'o, d{z) = -(32h/{2tl). Con- 
sider a two-step periodic dispersion map (see Fig. 1): 
d{z) — (d) + d{z), where d{z) = di for < z + < Li 
and d{z) = d2 for Li < z + nL < L, L = Li + L2 
is the dispersion map period, (d) is the path-averaged 
dispersion, di,d2 are the amplitudes of dispersion vari- 
ation subjected to condition diLi + ^2^2 = and n is 
an arbitrary integer number. Path averaged dispersion is 
assumed to be smaU, \{d) \ ^ \di\, 1^2 1. 

Optical pulse experiences strong oscillation as a func- 
tion of z on the period of dispersion map L. These os- 
cillations are caused primarily by linear effects which are 
described by the linear part of Eq. |j2Jl. 

A typical power of pulse is about a few mW in most 
optical fiber systems. As a result the nonlinearity in Eq. 
^ can be treated as a small perturbation on scales of 
typical dispersion map period L 50km because the 
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FIG. 1: Dispersion d{z) as function of propagation distance 
z. 



characteristic nonlinear length of the pulse is large: 
Zni ^ L, where Z^i — and p is the typical pulse 

amplitude. As nonlinearity is small perturbation on dis- 
tance L, it is convenient to introduce the new variable?/' : 
= ue"^'I£^/2di-')d^\ Here i>{uj) = jT^i^{t)e'^'dt is 
Fourier transform of if) and similar to u. ip is the slow 
function of z on the scale L because all fast linear de- 
pendence on z is included into fast varying exponent 



Fourier transform of Eq. Q gives: 



c{z) 



exp 



iA I d{z')dz' 

lLi/2 



Xl('{uJl)lp{uj2)'>P* i^3)H^l + ^2 ^ ^3 — i^)dLUidLU2dLU3,{3) 



Because ifj is the slow 



where A = cjj + — a;| — 
function one can approximate its derivative as '>p{z)z — 
[i('{z + L) — ip{z)]/ L and integrate Eq. ^ over dispersion 
map period L neglecting the dependence of ^ on z in 
the nonlinear term. This averaging procedure results in 
path-averaged Gabitov-Turitsyn ^ |^] model: 



iipzii^) - Lo'^{d)'il) 
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-iJj{uJi)iI){llI2) 



(27r)2 

X'0*(ljJ3)(5(cJl + LU2 — W3 — Lu)dujiduj2duj3, 

Co = (ciLi +C2i2)/i, (4) 

where s = diLi is the dispersion map strength, and 
ci, C2 are the values of c{z) for < z + nL < Li and 
Li < z + nL < L, respectively. It is set below, without 
loss of generality, that cq = 1 because one can always de- 
fine typical power Pq in such a way to insure cq — \. Note 
that Gaitov-Turitsyn Eq. can be also obtained by aver- 
aging of the Hamiltonian, H — J [d{z)\'4u\'^ — \i(^\'^ /2\dx, 



of NLS ^ over dispersion map period |^ . The numer- 
ical solutions of the Gabitov-Turitsyn model were com- 
pared with simulations of full NLS Q which showed good 
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agreement [ill Eoj . Returning to time domain in Eq. Q 
one gets |20l |: 

+ {d)4,tt - 2^ / Ci(^)^(^i + t)i'{t2 +t)x 

ip*iti+t2 + t)dhdt2 = 0, (5) 

where Ci(x) ~ ^^^^dx (note difference in definition 
of Ci{x) in comparison with |20j). 

Equation for DM soUton solution, V = A{t)e^^^ {A is 
real), of the Gabitov-Turitsyn Eq. ^ takes the following 
form: 

-\A+{d)Att = j r{tut2,t)dtidt2, 

r{ti,t2,t) = Ci{—)A{ti + t)A{t2 + + t2 + 0,(6) 

s 

where Ci{x) = J^cosx/xdx. 

Eq. 10 is used throughout this Article Eq. to study 
DM soliton. 

DM SOLITON TAILS 

The previous studies HSIllllSIllillllillillil 

illillillil have essentially focused on finding of DM 
soliton width, amplitude and soliton shape near soliton 
center. It was found that Gaussian ansatz, 



(a) 



Acauss = pexp ( - ^i^), 



(7) 



where p, [3 are real constants, is a rather good approxima- 
tion for the DM soliton solution near soliton center 
Solution of Eq. © by iterations allows to find soliton am- 
plitude and width analytically with accuracy ^1% (see 
Refs. 21, 27J). Note that parameters s, id), A uniquely 
determine the numerical form of DM soliton |27l| . E.g. for 
{d) =Q,f3s = 2.393 . . . , X=p'^x 0.482 . . . [STThus dis- 
persion map strangth, /3s, is fixed and not small for DM 
soliton. However an important question about asymp- 
totic behavior of soliton tail for large time t was only 
marginally addresses in the past ,2^i22i]. Soliton tails are 
responsible for interaction of pulses launched into optical 
fiber which essentially limits bit-rate capacity of optical 
line and makes finding DM soliton asymptotic behaviour 
an important practical and fundamental problem. DM 
soliton tails behaviour is the main subject of this Article. 

Solid curves in Fig. 2a and Fig. 2b show high preci- 
sion numerical solution of Eq. © for zero (Fig. 2a) and 
nonzero (Fig. 2b) values of averaged dispersion. Numer- 
ical technique for solving Eq. ® was developed in Refs. 
[27ll28j|. Analysis of numerical solutions in Fig. 2 allows 
to make hypothesis that at leading order the asymptotic 
behaviour of A is given by 




Time (t) 



(b) 




FIG. 2: DM soliton shape (curve 1) versus Eq. (4) (curve 2) 
for (a) (d) = 0, s = 1, A = 1. and (b) for (d) = 0.1, s = 
1, A = 1. A{t) is an even function. 



where oq, b are constants and /(i), a{t) are slow functions 
of t. First it will be shown that an ansatz ^ is consistent 
with Eq. © at leading order of \/\t\ provided ap = 
l/2s for f{t) — const, a{t) = const and later it will be 
found an expansion of f{t), a{t) over small parameter 
l/\t\. That result was announced without derivation in 
Ref. [i^l for particular case (d) = 0. Here general case 
which includes both (rf) = and (d) ^ is considered. 

First necessary condition for consistency of the ansatz 
© with Eq. © is that a substitution of the envelope 
of soliton solution tails, A^nv ^ exp(— 6|t|), into right 
hand side (r.h.s.) of Eq. © should recover the same 
exponential dependence on t after integration over ti,t2 
variables in Eq. ©. Asymptotic behaviour for — > oo 
of kernel in r.h.s. of Eq. © is given by: 



Ci{x) = sin{x)/x + 0{l/x^) 



(9) 



^{t) — f{t) cos{t'^[ao + a{t)]} exp{—b\t\), (8) and Ci(fif2/s) produces at leading order the contribution 
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FIG. 3: Integration plane (ti,t2) for Eq. (3). t > 0. Char- 
acteristic layer is located along AB inside the dotted curve. 
lushnikovfig3.eps 



to phase of A{t) only but not to envelope Aenv The 
contribution of other factors in r{t, ii, to A^nv is given 
by 

A{ti + t)A{t2 + t)A{ti +t2 + t) 
~ exp ( - b[\ti +t\ + |<2 +t\ + \ti +t2+ (10) 

and decays much faster than e"*"!*' outside the triangle 
ABC (shown in Fig. 3 for t > 0). Triangle ABC is 
defined in (ti,i2) plane and determined from condition 
h + t > 0, t2 + t > 0, tl + t2 + t < 0. However, in- 
side the triangle A{ti + t)A{t2 + t)A{ti + 12 + 1) ~ e'''^*^ , 
which recovers the exponential decay as it appears at 
left hand side (l.h.s.) of Eq. ©. One can conclude that 
Aenv ~ e^''l*l is compatible with Eq. © and integra- 
tion in plane (^1,^2) inside an area outlined by dashed 
curve in Fig. 3 gives leading order contribution to r.h.s. 
of Eq. © for \t\ 00. The area inside the dashed 
curve includes triangle ABC and characteristic layers of 
the width 0(t°) which surround that triangle. Integra- 
tion over the area outside of dashed curve in Fig. 3 
gives exponentially small contribution of r.h.s. of Eq. 
®. E.g. A{h + t)A{t2 + t)A{ti i2 + t) =i e-sf-l*! for 
ti~i/2, t2 ^ t/2. This result is immediately confirmed 
by numerical calculation of integral in r.h.s of Eq. © 
over the area inside the dashed curve for DM soliton so- 
lution and comparing result of integration with numerical 
integral over the whole {ti, 12) plane. Numerical integra- 
tion was performed using numerical procedure of Ref. 
[27I 1. Note the important fact that exponential envelope, 
Agnv , is compatible with Eq. ^ for any positive value of 
b so that specific value of b, corresponding to DM soliton 



can be only obtained from more detail analysis of Eq. 

Second necessary condition is that the substitution of 
fast oscillating function Aosc ^ cos {P[ao + a{t)]) into 
r.h.s. of Eq. © should allow to reproduce the same 
oscillations as in l.h.s. of Eq. for |<| 00. Oscillating 
part of integrand r(i, tl, i2) of r.h.s. ofEq. Q consists of 
product of one sin function (from asymptotic of oscillat- 
ing kernel ©) and three cos functions. Representing sin 
and cos through imaginary exponents one get a sum of 16 
terms with different values of argument. 14 of them are 
always oscillate as a quadratic function of ti , t2 and thus 
give vanishing contribution after integration over ti, t2- 
Rest two terms can be made nonoscillating by appropri- 
ate choice of hq. Thus integrand r(t, ii, t2) of r.h.s. ofEq. 
lO at leading order, neglecting a{t) inside the triangle 

ABC, can be represented as r{t,ti,t2) ^ exp {itit2/s + 

iaolih+tf + {t2 + tY - {ti+t2 + tf]-b\t\^ -c.c. + f.o. = 

e:xjp{itit2/s — 2iaQtit2 + iaot^ — b\t\) — c.c. + f.o., where c.c. 
mean complex conjugation and f.o. designates 14 fast os- 
cillating terms. Oscillations (as a quadratic function of 
tl, t2) of argument of two rest terms vanish identically 
provided 

ao = l/(2s), (11) 

which is in excellent agreement with numerical solution 
of DM soliton ©. Next step is to take into account 
nonzero value of a{t) in Eq. JHjl. One can series expand 

a{t)^ai/\t\+a2/t^ + 0{l/\t\^) (12) 

and find that term ai|t| in argument of cos of Eq. 
causes oscillations of r(t, ti,t2) as a linear function of 
ti,t2. It is impossible to get rid of these oscillations for 
any value of ai so one can conclude that contribution 
of integration of r(t,ti,t2) inside the triangle ABC van- 
ishes meaning that leading order contribution to r.h.s 
of Eq.jni can only come from characteristic layers of 
the width 0(t°) which surround faces AB, BC and AC 
of triangle ABC. In the neighborhood of each of these 
faces either ti + t (face AC) or t2 + t (face BC) or 
^1+^2 + ^ (face AB) is ~ t° so asympotic expression 
is not valid but rather exact solution of Eq. © 
should be used for respective factors A{ti + t), ^(t2 + t) 
or A{ti + t2 + t) in r(t, ti,t2). Still there are fast os- 
cillations of r{t,ti,t2) as a function of ti,t2 for inte- 
gration inside the characteristic layers in directions par- 
allel to faces BC and AC: A{ti + t)A{ti + t2 + t) ^ 
gmi[2t+2ti+t2]-6|t2| ^ ^. ^ J Q (along BC) and ^(t2 + 

t)A{ti +t2+t) ^ giai[2t+2ti+t2]-b|ti| ^_ ^ ^ J ^ (j^loj^g 

AC). Thus integrations in these 2 layers give a vanish- 
ing contribution and only an integration in characteristic 
layer in the neighborhood of face AB gives leading order 
contribution. To check that result the numerical integra- 
tion of integrand r{t,ti,t2) with super Gaussian modi- 
fication factor: f(t,ti,t2) = r(t, ti, t2)e-(l*+*i+*2l/*o)'''"' , 
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n 1, 1 <^ to t was performed and it was found that 
J fit,h,t2)dtdt2 = Jr{t,h,t2)dtdt2 + 0{l/\t\) which 
confirms that only integration inside the characteristic 
layer around AB (that layer is schematically shown in 
Fig. 3 as the area inside the dashed curve) gives leading 
order contribution. 

Analysis of an integration inside the characteristic 
layer AB is convenient to be performed in new variables: 



^ = t + h+t2, 

V^h- t2, 



(13) 
(14) 



where axes ^ and r] are perpendicular and parallel to 
face AB, respectively (see Fig. 3). Face AB is located 
exactly at axis ^ = 0. Keeping in mind that the width 
of characteristic layer is ~ i° ~ ^ while its length is 
^/2\t\ ~ ry one gets that the ratio of typical values of ^ 
and ?7 is - l/\t\ and h + t = ^±2 + 0(^0)^ t2 + t = 
^ + 0(^0), = - *^ + 0(t"), t2 = + 0(i°). So 
that Eq. @, using Eq. lO, is reduced at leading order 
of l/\t\ to 



Xgiid)) + {d)ty.s^] fit) cos + a,\t\ + 03) 



r 

X exp{-b\t\) ^ J ^ J drj 



15) 



(r,-i)(77 + i) 


X sin ([e^ + t^]/2s + ai[C+ \t\] + 2a2) exp {~b( - b\t\), 

where the function g(x) is used, g{0) = 1 and g{x) = 1 
for X 0, because term with A in l.h.s. of Eq. (|15|l is 
of leading order only if (d) = 0. There is a separation of 
the integration variables ^, 77 in r.h.s of Eq. and one can 
perform integration over ^ to get reduced equation: 



Xgiid)) + {d)tys^ /(i)-Ci 



v)fiv) 



it - 77)77 



(16) 



where 77 — it is assumed without loss of generality 
that t > 0; the real constant Ci is given by 

Ci= f ^^(0 exp he /2s + + 02) - 



Ci=i?e(Ci), (17) 

and Ci — i?e(Ci) is the compatibility condition for Eq. 
(|15|l meaning that phase of cos in both r.h.s. and l.h.s. 
(after integration over ^) is the same and equals to 02 
for |t| = 0. It is important that exact (not asymptotic) 
expression for Ait) should be used in evaluating integral 
in Eq. 1)171) which means that asymptotic tails © are 
determined by strongly nonlocal effects because all time 
scales contribute to Ci . 

The integral term in Eq. (|16|) is the convolution inte- 
gral and using Laplace transform, fp = dt/(t)e~^'. 



for function fit) = /(t)/t one gets an ordinary differen- 
tial equation: 



Aff((rf)) 



dfp 
dp 



id) 



1 d^fp 
dp^ 



-Cifp, 



(18) 



which has the following solution: fp ~ giid)) Ci{p+C2) 
SOTT-TTTTT-TTr, whcrc C2 is an arbitrary real constant. 

Laplace transform results in: fit) = 



'Ciip+C^r ' 
Inversion of 



tgiid))\ + ?>Qt^id) e~^2t/Ci. For arbitrary sign of t one 

should replace thy \t\ in this expression. Factor e^*-^^'*' 
gives renormalization of constant b in Eq. (jSJ, which re- 
flects a freedom of choice of b in Eq. © so C2 could be 
set, without loss of generality, to zero and expression for 
fit) takes the following final form: 



A|t| 



QO \ 1 + 13 

fit) = '-^ for id) = 0, f(t) = — ^ for id) + 0. 

(19) 

Eqs. (|12)1 . ()19(l have four real constants Ci, 6, a\ and 
a2 which can not be determined from leading order ex- 
pansion developed in this Article because all time scales 
contribute to numerical values of these constants. One 
can easily find constants Ci, 6, Oi, 02 from a fit with nu- 
merical solution of Eq. ©. Phase of cos in Eq. JSJ) is 
determined from zeros of the numerical solution of A(t) 
keeping in mind that coefficient oq = l/2s is known an- 
alytically. After the constants ai, 02 are found, one can 
determine C\ and b numerically from a two-parameter fit 
with envelope of DM soliton tail oscillations. Fig. 2a, b 
show a good agreement between the numerical solution 
(curve 1) and Aasymp^f)'^ dependence (curve 2). Fig. 2a 
corresponds to s = 1, A = 1, (d) = and numerical val- 
ues Ci = 0.08357, b = 3.0452, ai = 1.4136, 02 = 1.5102. 
Fig. 2b corresponds to s = 1, A = 1, (c?) =0.1 and numer- 
ical values Ci = 13.818, b = 3.2483, oi = 0.1443, 02 = 
—0.1276. Ci, 6, oi, 02 are obtained from the above de- 
scribed fitting procedure. A similar plot for zero average 
dispersion (d) — was given in Ref. [l^l in Fig. 3. Note 
however that curve 2 in Fig. 2b converges to curve 1 for 
larger t compare to the case (d) = in Fig. 2a which re- 
flects the fact that for intermediate values of time, i ~ ii, 
(d)tf = A, there is a transition from the solution of f{t) 
in Eq. (O for (d) = to the solution for (d) 9^ 0. If 
(d) is large enough so that ti ^ DM soliton width, then 
there is a transition from oscillating tails of DM soliton 
to nonoscillating exponential tails of NLS with constant 
dispersion (d). NLS soliton width, (A/(d))-^/^, is much 
larger than DM soliton width for small (d) . 



CONCLUSION 

In summary, the asymptotic, \t\ 00, of DM soli- 
ton tails is given by Eqs. © , IIIll , , (CHI) which is 
in agreement with numerical solution of Eq. ©. The 
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pre-exponential factor H19|l is qualitatively different for 
zero and nonzero path-averaged dispersion. An exponen- 
tial decay of DM soliton envelope, which is much slower 
compare to Gaussian-type decay near soliton center, sug- 
gests that interaction of different DM solitons is strong 
and could affect bit-rate capacity of optical lines. A de- 
tailed investigation interaction between DM solitons is 
outside the scope of this Article. Another possible direc- 
tion of future research is to study the effect of losses on 
DM soliton tails as well as DM soliton with nonlinearity 
compensation Qualitatively oscillations are still of 

the form (jSJ but further research is needed. 
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